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ON  THE  STABILITY  OF  VORTEX  MOTIONS 
IN  THE  PRESENCE  OF  MAGNETIC  FIELDS 

INTRODUCTION 

In  an  important  paper  on  the  hydrodynamic  and  hydromagnetic  stability  of  nondissipitive  swirling 
flows,  Howard  &  Gupta  (1962)  presented  many  interesting  aspects  on  stability  characteristics  of  vortex 
motion  with  and  without  axial  velocity.  Even  though  they  restricted  themselves  mostly  to  homogene¬ 
ous  fluids  subject  to  axisymmetric  perturbations,  their  generalization  of  the  Richardson  criterion  and  of 
the  semi-circle  theorem  to  axisymmetric  steady  flows  presented  a  relatively  simple  insight  into  the  type 
of  parallel  flows  complicated  by  the  cylindrical  geometry.  Acheson  (1972,  1973)  later  extended  the  ana¬ 
lyses  into  the  same  type  of  flows  subject  to  non-axisymmetric  perturbations.  A  review  paper  on  this 
subject  was  given  by  Acheson  St  Hide  (1973). 

By  considering  a  homogeneous  fluid  rotating  uniformly  in  a  radius-dependent  magnetic  field, 
Acheson  (1972)  derived  some  sufficient  conditions  for  stability  of  the  flow.  The  difference  between  the 
axisymmetric  and  nonaxisymmetric  modes  were  brought  out  by  assuming  the  wavelengths  in  the  radial 
direction  to  be  small  compared  with  the  radius.  In  addition,  the  stability  phenomenon  that  all  non¬ 
axisymmetric  unstable  waves  must  propagate  against  the  basic  rotation,  i.e.,  the  westward  drift,  was 
proved  to  prevail  with  a  restriction  on  the  axial  and  azimuthal  components  of  the  magnetic  field.  With 
the  help  of  the  Boussinesq  approximation,  the  author  later  (Acheson  1973)  generalized  the  westward 
drift  phenomenon  to  heterogeneous  fluids  rotating  differentially.  Also  based  on  this  westward  drift,  a 
quadrant  theorem  reminiscent  of  the  semi-circle  theorem  encountered  in  two-dimensional  stratified 
flows  was  derived  for  slow  amplifying  waves. 

It  should  be  emphasized,  as  also  noted  explicitly  by  Acheson  (1973),  that  the  results  obtained  in 


his  paper  are  restricted  to  perturbations  with  non-zero  axial  wave  numbers.  Any  attempt  to  infer  that 
Manuscript  approved  October  12,  1982. 


the  results  can  be  applied  to  perturbations  with  zero  axial  wave  numbers  (azimuthal  modes)  may  lead 
to  incorrect  conclusions  on  both  the  hydrodynamic  and  hydromagnetic  problems.  This  argument  will 
be  demonstrated  by  examining  the  azimuthal  modes  in  this  paper. 

In  the  present  investigation,  stability  analyses  are  performed  on  a  general  type  of  vortex  flow  with 
varying  density  in  the  presence  of  an  axial  and  an  azimuthal  magnetic  field.  A  sufficiency  condition  for 
stability  is  derived  and  compared  with  the  results  obtained  from  other  methods.  The  dual  role  played 
by  the  rotating  velocity  on  the  stability  of  vortex  motions  is  separated  and  revealed  by  proper  transfor¬ 
mations.  While  the  shear  effect  conveyed  by  the  velocity  gradient  at  the  shear  layer  always  destabilizes 
the  flow,  the  centrifugal  force  generated  by  the  fluid  rotation  stabilizes  or  destabilizes  the  flow  depend¬ 
ing  on  whether  or  not  the  density  increases  radially  outwards. 

For  vortex  motions  subject  to  azimuthal  disturbances,  a  magnetic  field  always  has  a  stabilizing 
effect  regardless  of  its  detailed  distribution.  As  an  extension  of  the  semi-ellipse  theorem  in  stratified 
rotating  flows  (Fung  1982),  an  upper  bound  on  azimuthal  amplifying  waves  is  derived  with  a  restriction. 
The  restriction  suggests  that  the  semi-circle  in  the  complex  velocity  plane  does  not  provide  a  bound  on 
all  the  growing  azimuthal  modes.  Furthermore,  the  upper  bound  demonstrates  that  unstable  waves  do 
not  necessarily  drift  westward  for  the  hydromagnetic  or  non-hydromagnetic  case.  For  uniformly  rotat¬ 
ing  flows  with  constant  angular  Alffen  velocities,  all  azimuthal  unstable  waves  must  lie  on  a  semi-circle 
independent  of  the  density  distributions.  Exact  solutions  for  some  special  flow  profiles  are  obtained  to 
support  the  derived  stability  criteria. 

MATHEMATICAL  DERIVATION 

Consider  a  vortex  flow  with  a  velocity  V  and  a  magnetic  field  H  to  be  confined  within  the  annular 
region  (r,  9,  z)  between  two  rigid,  infinite  and  coaxial  cylinders.  The  fluid  having  an  inhomogeneous 
density  p*  is  assumed  to  be  inviscid,  incompressible,  and  non-heat-conducting.  When  gravitational 
forces  and  dissipation  effects  due  to  viscosity,  magnetic  resistivity,  and  thermal  diffusivity  are  neglected, 
the  governing  equations  for  the  flow  are 
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p*  -TT-  VQ  +  -f-  (ff  V)  H  (1) 

Dt  w 

V  •  ff  -  0  (2) 

M-Vx(ffxff)  (3) 

dt 

V  •  ff  -  0  (4) 

^-0  (5) 

Dt 

where  p  denotes  the  magnetic  permeability.  The  total  pressure  Q  (including  the  magnetic  pressure)  is 
related  to  the  hydrodynamic  pressure  P  as  follows: 

Q  -  P  +  Iffl2.  (6) 

The  boundary  conditions  for  the  system  governed  by  Eqs.  (1)  to  (5)  are  those  of  perfectly  conducting 
rigid  walls. 

The  flow  to  be  considered  has  a  steady-state,  radius-dependent  profile  as  follows:  (1  (r)  is  the 
angular  velocity,  tV(r)  the  axial  velocity,  Ht(r)  the  azimuthal  magnetic  field,  Hz(r)  the  axial  magnetic 
field,  Q0(r)  the  total  pressure,  and  p0(r)  the  density.  Let  the  flow  be  perturbed  as  follows: 

ff  -  ff  [u,  r(l(r)  +  v,  W(r)  +  w), 
ff  -  ff  [ft,,  H,(r )  +  it „  H.(r)  +  it,], 

Q  -  Qo(r)  +  q, 

P*“Po(r)+p.  (7) 

We  further  introduce  the  periodic  solutions 

4  —  ^(r)  Exp  U(kz  +  m9  —  «f)l  (8) 

such  that  the  azimuthal  wave  number  m  is  an  integer,  the  axial  wave  number  k  is  real  and  positive,  and 

the  circular  frequency  at  -  »,  +  lot,  is  complex.  Within  the  framework  of  the  normal  mode  method, 

the  linearized  equations  for  the  flow  described  by  Eqs.  (1)  to  (S),  subject  to  small  perturbations,  are 

given  as  follows: 
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(12) 

(13) 

(14) 

(15) 

(16) 
(17) 


H  d 

where  N  ~  kW  +  mil  —  u  is  the  Doppler-shifted  frequency,  N„  —  kH,  +  m  — -,  D  “  — ,  D*  —  Z) 

r  dir 

+  -y  and  Z).  —  Z>  -  H  The  characteristic  that  Eqs.  (3)  and  (4)  represent  only  three  independent  par¬ 


tial  differential  equations,  under  the  present  assumption,  is  reflected  in  Eqs.  (13)  to  (16). 


If  we  define  the  angular  and  axial  Alfv6n  velocities  as 


Ha 


and  the  Alfven  frequency  as 


NA  •  *  WA  +  «  &A’ 

Eqs.  (9)  to  (17)  can  be  combined  into  two  independent  differential  equations  as  follows: 


where  the  Rayleigh-Synge  discriminant  and  the  Alfv6n  discriminant  are  respectively  defined  as 
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The  equation  governing  stability  of  the  flow  can  be  obtained  by  combining  Eqs.  (18)  and  (19)  as  fol¬ 
lows: 


D  I p0E(N2  -  A rj)Z>**]  -  \D. 


2mpoE 


-  NAClA) 


+  pol AT2  -  Nj  -  4>  +  VA  +  4H2  -  4 t^E 


(AW  -  NAClA)2 


N2  -  N} 


U  -  o 


where 


and 


V  N 


(20) 


(21) 


E  - 


m2  +  fcV ' 

The  boundary  conditions  for  the  above  equation  are 


*(r,)  -  *(r2)  -  0  (22) 

where  r,  and  r2  are  locations  of  the  solid  boundaries.  If  discontinuities  in  the  flow  profile  exist  within 

the  flow  domain,  the  matching  conditions  obtained  by  integrating  Eqs.  (18)  and  (19)  across  the  discon¬ 
tinuity  surface  are 


<-£)-0 
N 

<«>-  <P0r( n2 -  n2)>  -  o 


(23) 

(24) 


where  <$>  —  <t>(n+o)  ~  4hr- o>  represents  the  jump  condition  across  the  discontinuity  surface  at 
r-R. 


STABILITY  CRITERIA 


We  will  investigate  the  general  stability  characteristics  described  by  Eq.  (20).  Due  to  the  com¬ 
plexity  of  the  presence  of  all  the  four  components  in  the  velocity  and  magnetic  fields,  general  stability 
criteria  only  in  terms  of  the  flow  profiles  are  difficult  to  derive.  However,  stability  conditions  for  some 
particular  cases  can  still  be  observed  and  discussed  in  the  following. 
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(1)  Snttcfent  Condition  for  Stability 


Let 
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(25) 


Equation  (20)  is  then  transformed  to 
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Multiplying  the  above  equation  by  r<ft,  where  the  quantity  with  a  bar  stands  for  the  complex  conjugate, 

and  integrating  the  resultant  equation  over  the  flow  domain,  we  obtain 


J*Po|(V - (E\D*$\2  +  \<t>\ 2)rdr 

(l  -  +  TT 

fll  <,£l DN)! 

J  (4  N 


-  —  •  -4il2  +  4*2£ 


N2  j  rN2 
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with  its  imaginary  part  equal  to 
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Since  the  first  three  integrals  are  all  positive  definite,  no  solution  corresponding  to  m,  &  0  exists  if  the 
sum  of  the  quantities  in  the  last  integral  in  Eq.  v28)  is  positive.  Therefore  we  can  conclude  that  the 
flow  is  stable  if 


<P-VA-  4(ft2  -I-  0^)  +  (kWA  +  mClA)2 


r2(kDW+mD(l)  2m(l  kDW  +  mDCl 

\  ,  (kWA  +  m(lA)2 

m 2  +  kV  r  4 

1  kW  +  mSl  —  &>  | 2 

anywhere  within  the  flow  domain.  Equation  (29)  can  be  rearranged  to 


0  (29) 


(1  +  s2){4>  -  VA  -  4H  j  +  tfj)  -  4n2  1 


[DW  +  s(rD(l  +  4ft)]2  + 


( DW  +  srDCl ) 


Na_ 

I  AM 


(30) 


where  s  —  m/  kr.  It  can  be  seen  immediately  that  Eq.  (30)  will  be  violated  when 


<D  -  VA  -40^  +  N}  <  4n2/(l  +  s2).  (31) 

Equation  (30)  is  reminiscent  of  the  stability  criterion  obtained  by  Fung  &  Kurzweg  (1975)  in  their 

study  on  the  stability  of  swirling  flows  with  radius-dependent  density.  Like  the  role  played  by  the 

Rayleigh-Synge  criterion  (a  condition  for  centrifugal  stability)  in  their  stability  criterion,  Eq.  (31)  can  be 

viewed  as  a  precondition  to  the  sufficiency  criterion  in  Eq.  (30).  However,  this  sufficiency  criterion  for 

zero  magnetic  forces  does  not  reduce  to  the  one  obtained  by  Fung  &  Kurzweg  (1975).  This  is  a  well- 

known  paradox  stemming  from  the  fact  that  the  lines  of  magnetic  flux  for  a  fluid  with  zero  resistivity 

are  permanently  attached  to  the  fluid  (Chandrasekhar  1961).  The  attachment  is  contained  in  the 

second  and  the  third  integral  of  Eq.  (28).  If  we  assume 


Na  -  kWA  4-  mClA  -  0 

throughout  the  flow  domain,  Eq.  (28)  can  be  written  as 


Po(£ \D*$\2  +  l^l2)  rdr 


(32) 


{/ - 

+  /  |/2  +  /j  +  *  -  vA  -  4(n2  +  n \)  -  _  UMlL  J  «  0. 


The  positive  definite  integrands 


h- 


4m2flj 

m 2  +  kV 


(33) 


(34) 
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and 


(35) 


m2  +  frV 

which  were  originally  not  admitted  in  the  sufficiency  condition  obtained  from  the  second  and  third 
integral  in  Eq.  (28),  are  now  recovered.  The  resultant  condition  for  stability  now  reads 


(1  +  s2)  (<t>  -  VJ 


>4. 


[DW  +  s(rD(l  +  40)]2  4'  (36) 

Note  that  Eq.  (36)  is  valid  only  with  the  constraint  described  by  Eq.  (32).  It  is  remarkable  that  even 

though  the  constraint  involves  the  magnetic  forces  in  both  the  axial  and  azimuthal  direction,  the 

sufficiency  condition  described  in  Eq.  (36)  does  not  depend  on  the  axial  magnetic  forces.  Furthermore, 

if  we  compare  the  numerator  in  Eq.  (30)  with  that  in  Eq.  (36),  a  destabilizing  term 

-4[fl2  +  (1  +  s2)fl  21  appearing  in  the  former  disappears  in  the  latter  as  NA  approaches  zero.  Though 

this  paradox  suggests  a  destabilization  of  the  flow  by  the  magnetic  forces,  one  should  keep  in  mind  that 

both  the  sufficiency  conditions  in  Eqs.  (30)  and  (36)  only  represent  a  bound  on  stability  and  are  by  no 

means  to  give  the  final  stability  conditions. 


Since  the  sufficiency  condition  for  stability  obtained  in  Eq.  (29)  or  (30)  contains  the  unknown  fre¬ 
quency  <o,  a  more  general  view  of  flow  characteristics  in  terms  of  given  flow  profiles  is  therefore 
difficult  to  see.  In  the  following  we  will  further  examine  the  flow  characteristics  using  other  methods 
and  compare  the  results  with  the  derived  sufficiency  conditions  for  some  particular  cases. 

(2)  Comparisons  with  Different  Modes 

Three  different  types  of  perturbations  are  to  be  examined  and  compared  with  the  earlier  obtained 
sufficiency  condition. 


A.  Axisymmetric  modes  (m  -  0) 

For  simplicity,  we  will  neglect  the  axial  magnetic  flux,  i.e.,  WA  «■  0  throughout  the  flow  domain. 
Equation  (20)  under  the  present  conditions  reads 

D  \p0(W-  ck)2D**]  -po  l k2(W-  c*)2-4>  +  *4]*-0  (37) 
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where 


co  .  . 

cfc  -  —  -  cr*,  +  icu 

is  the  axial  phase  velocity.  The  constraint  in  Eq.  (32)  is  now  satisfied  and  the  loss  of  the  stabilizing 
effect  by  the  presence  of  the  axial  magnetic  flux  is  recovered.  The  sufficient  condition  for  stability 
reduced  from  Eq.  (36)  becomes 


4  >  -r 


(38) 


where 


<p-V, 


(39) 


*  ( DW )2 

which  is  reminiscent  of  the  Richardson  number  in  two-dimensional  stratified  flows.  Condition  (38)  was 
first  derived  by  Howard  &  Gupta  (1962)  for  homogeneous  fluids.  In  the  absence  of  axial  flows,  Eq. 
(37)  together  with  the  boundary  conditions  forms  a  Sturm  Liouville  system  for  which 


0  (40) 

is  a  necessary  and  sufficient  condition  for  stability.  The  above  condition  was  first  derived  by  Michael 

(19S4)  for  homogeneous  fluids  and  will  be  called  the  generalized  Michael  condition.  An  alternative 

way  to  examine  the  influence  of  flow  quantities  on  stability  can  be  obtained  as  follows. 


Multiplying  Eq.  (37)  by  nfi  and  integrating  the  resulting  equation  throughout  the  flow  domain,  we 
obtain 

/  (W  —  ck)2  Xkdr  -  /  p0(«&  -  V4)  1*1 2rdr  -  0  (41) 

where 


X*  -  p0  (|Z)*ifr|J  +  fc2|i|r|2)r  >  0- 

Solving  Eq.  (38)  for  ck  results  in 


where 


J  WXkdr  ±  JL  Xkdr  fpo  (»-  ^)  1  *l>\2rdr  —  6 


(42) 


«*  -  /  X*dr/  W*Xkdr  -  (/  WXkdr)2  >  0 
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for  all  values  of  W  resulting  from  the  Schwarz  inequality.  It  follows  from  Eq.  (42)  that  instability  of 
the  flow  is  expected  when 

—  8*  +  J  Xkdrf  Po  (4>  —  ¥4)  \\ltVrdr  <  0.  (43) 

It  is  obvious  from  the  above  equation  that  the  axial  velocity  always  destabilizes  the  flow  except  for  con¬ 
stant  W  where  8*  =  0.  The  second  term  in  Eq.  (43)  represents  the  contribution  from  the  Rayleigh- 
Synge  and  the  Alfven  discriminants.  Instability  will  occur  when  the  generalized  Michael  condition  is 
violated.  It  is  interesting  to  note  that  even  though  both  Eqs.  (38)  and  (43)  involve  similar  arguments 
on  the  axial  velocity  and  on  the  generalized  Michael  condition,  they  represent  different  bounds  on  flow 
stability.  It  is  obvious  from  Eq.  (43)  that  violating  the  generalized  Michael  condition  automatically 
leads  to  instability  of  the  flow.  However,  such  a  conclusion  can  not  be  drawn  directly  from  Eq.  (38) 
since  violating  the  sufficiency  condition  does  not  necessary  lead  to  instability. 


B.  Azimuthal  modes  ( k  —  0) 


The  sufficiency  condition  for  k  -  0  reduced  from  Eq.  (29)  reads 

_ ^ 

l+fii/(ft-cm)J  "  4 

r(l2Pp0  —  rP(p0Sl%)  +  (m2  -  4)p0fl  j 

m"  PoirPd)2 

where 


(44) 

(45) 


cm  -  “  -  <W  +  icmi 

m 

and  is  the  angular  phase  velocity.  An  analogy  can  be  drawn  between  the  two  sufficiency  conditions  in 
Eqs.  (38)  and  (44)  except  in  the  latter  the  shear  effect  is  produced  by  the  angular  velocity  rather  than 
the  axial  velocity.  Equation  (44)  suggests  that  the  angular  velocity  plays  a  dual  role  in  flow  stability: 
the  angular  velocity  itself  stabilizes  the  flow  while  its  gradient  destabilizes  the  flow.  This  characteristic 
can  further  be  observed  in  the  following  analysis.  The  influence  of  the  azimuthal  magnetic  field  on 
flow  stability  is  difficult  to  see  because  the  sufficiency  condition  involves  the  complex  phase  velocity. 
An  alternative  way  to  examine  the  flow  characteristics  is  given  as  follows. 
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Equation  (20)  for  k  —  0  is  written  as 
Dlpo^iN1  -  m2Ci\)Dm>b\ 

-  {2wi>lp0(/Vft  -  mfi^)]  +  p0m2( N2  -  4>  +  4fl2  +  -  m2ftj)}0  -  0.  (46) 

Here 

N  ~  mCl  —  a. 

Multiplying  Eq.  (46)  by  n)i  and  integrating  the  resultant  equation  yield 

f  (N1  -  m2(l\)  Xmdr  +  f  { 2mWV£>(p0ft )  -  m2  [ ril2(Dp0 )  4-  rD(p^£l\)])r  \th\2dr  -  0  (47) 

where 

Xm  “  Pti(r2  |Z>*0l2  +  m2  |i|»|2)r  >  0.  (48) 

Further  we  utilize  the  transform 

Y*  “  Po ir2\D$\2  +  ( m 2  -  1)  |i|t!2lr  >  0  (49) 

Eq.  (47)  can  be  reduced  to  a  simple  quadratic  form  of  the  complex  azimuthal  phase  velocity  cm 

cif  Xmdr  -  2  cmJ  ft  Y  mdr  +  /  (ft2  -  Cl2A)Ymdr  -  0.  (50) 

Solving  (50)  for  cm  leads  us  to  a  conclusion  that  the  flow  will  be  unstable  (corresponding  to  c,  *  0) 

when 

-»«  +  /  UW i3  l#l 2drf  ft2Y mdr  +  /  Xmdr  f  Slfr mdr  <  0  (51) 

where 

S„  -  /  Ymdr  J  a2Ymdr  -  (/  SlYmdr)2  >  0 

for  all  values  of  ft  resulting  from  the  Schwarz  inequality.  The  first  term  of  Eq.  (51)  represents  the 
shear  effect  produced  by  the  angular  velocity  gradient.  It  is  obvious  that  this  term  always  destabilizes 
the  flow  except  for  uniform  rotation  where  8m  —  0.  Instability  in  this  case  is  of  centrifugal  origin 
because  of  the  absence  of  shear  layers  within  the  flow  domain.  The  second  term  (51)  represents  the 
density  variation  within  the  centrifugal  force  field  produced  by  the  rotation  of  fluid.  Positive  density 
gradients  stabilize  the  flow.  In  the  absence  of  magnetic  fields,  instability  automatically  occurs  when  the 
density  decreases  radially  outwards  regardless  of  the  detailed  profile  of  the  angular  velocity.  Even 
though  it  is  also  suggested  by  Eq.  (44)  that  negative  density  gradients  may  lead  to  the  violation  of  the 
sufficient  condition  for  stability,  instability,  however,  can  not  be  concluded.  Violating  Eq.  (44)  does 


II 


not  necessarily  lead  to  instability.  The  last  term  in  Eq.  (51)  represents  the  magnetic  influence  on  flow 
stability.  The  presence  of  the  azimuthal  magnetic  force,  regardless  of  its  detailed  distribution,  always 
stabilizes  the  azimuthal  disturbances.  This  argument,  which  cannot  be  seen  clearly  from  criterion  (44), 
wii'  '  supported  by  an  analytical  solution  to  the  governing  stability  equation.  Another  conclusion  we 
can  draw  from  Eq.  (51)  is  that  for  uniform  rotation  flow  stability  can  be  guaranteed  if  the  density 
increases  radially  outwards.  This  characteristic  is  also  independent  of  the  profile  of  the  magnetic  field 
and  cannot  be  observed  from  Eq.  (44). 


C.  Arbitrary  modes 


Because  of  the  complexity  of  the  arbitrary  perturbations,  we  will  ignore  the  influence  from  the 
magnetic  field  and  concentrate  on  the  contribution  from  the  density  and  velocities  of  the  flow.  Under 
these  conditions,  the  constraint  in  Eq.  (32)  is  satisfied  and,  without  losing  the  the  stabilizing  effect,  the 
sufficient  condition  for  stability  reduced  from  Eq.  (36)  reads 


(1  +  s2)  d> 


>1 


(52) 


[DW  +  sirDSl  +4ft)l2  4‘ 

This  condition  was  first  derived  by  Fung  &  Kurzweg  (1975)  in  their  study  on  heterogeneous  swirling 

flows.  Readers  are  referred  to  their  detailed  discussions  on  the  condition.  Reminiscent  of  the  statically 

stable  condition  for  density  encountered  in  two-dimensional  shear  flows,  the  Rayleigh-Synge  criterion 

(4>  >  0)  acts  as  a  condition  for  centrifugal  stabilities. 


In  the  following,  we  will  further  investigate  the  stability  characteristics  by  adopting  the  method 
used  in  the  axisymmetric  and  the  azimuthal  cases.  The  integral  equation  thus  obtained  is 

/  ***  -*/-•»  *-/(*  -  ^rj  -  »  o» 

where 

X  “  Po|  m2  U>**l2  +  I*  lJ)  r  >  0.  (54) 
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We  further  utilize 


„  p0 lr2\D<l>\2  +  (k2r2  +  m2  —  Dl^Plr  ^  n 

Y  ■■  _  r  ,  ^  U. 

m2  +  ik2/2 

Equation  (53)  is  transformed  to  a  simple  quadratic  form  of  a 

a2f X<fr  -  2a  f  (kfVX  +  mClY)dr  +  f  (k2W2X  +  IkWmClY  +  m2Sl2Y)dr 

Instabilities  are  expected  whenever 

-8 o  -  8,  -  8,  +  J*  (£>p0)  m?+l2r2  drf 


+  fXdrf  Po*  -2^72  l*|2rdr  <  0 


Here 


(55) 


(56) 


(57) 


5*  -  /x</r  J*  ifc^Xdr  -  (/  ikJKXrfr)2 
Sj,  -  f  Ydrf  m2tl2Ydr  -  (/  mCiYdr)2 
8,  -  2(J  Xdrf  kWmtlYdr  -  f  kfTXdrf  mSlYdr) 

8/-2[jVn'Y<*/  -/x*/ 

Since  X  and  Y  are  positive  definite,  it  follows  from  the  Schwarz  inequality  that  both  S^  and  Sn  are 
positive  definite  for  all  values  of  W  and  SI ,  representing  the  shear  effects  conveyed  by  the  gradients  of 
the  axial  and  tangential  velocities.  As  also  implied  by  condition  (52)  they  both  destabilize  the  flow 
except  for  constant  axial  and  angular  velocities  where  8  ^  «  8n  —  0.  The  third  and  fourth  terms  are 
respectively  the  influence  of  the  perturbation  directions  on  the  shear  effect  produced  by  the  velocities 
and  on  the  centrifugal  force  created  by  the  fluid  rotation.  Both  can  be  either  positive  or  negative, 
implying  that  they  can  either  stabilize  or  destabilize  the  flow.  This  stabilization  or  destabilization  will 
not  be  seen  until  solutions  to  the  governing  stability  equation  are  obtained.  The  fifth  term  in  (57)  is 
the  effect  of  density  variations  in  the  centrifugal  force  field.  Positive  density  gradients  stabilize  the 
flow.  The  last  term  in  (57)  involves  the  Rayleigh-Synge  discriminant  and  stabilizes  the  flow  if  the 
Rayleigh-Synge  criterion  (<b  >  0)  is  satisfied.  Condition  (57)  can  reduce  to  condition  (43)  for  the 
axisymmetric  case  and  to  condition  (51)  for  the  azimuthal  case  if  the  magnetic  fields  are  deleted. 
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(3)  Bounds  on  Unstable  Wares 


For  instability,  a  bound  on  the  growth  rate  can  readily  be  obtained  from  Eq.  (28).  Since  the  first 
three  integrals  in  the  equation  are  all  positive  definite,  it  follows,  as  an  opposite  to  the  sufficiency  con¬ 
dition  for  stability  in  Eq.  (29),  that  the  integrand  in  the  last  integral  must  be  negative  somewhere 
within  the  flow  domain.  This  leads  to  a  bound  on  the  growth  rate  such  that 

2  ( _ N}(DW  +  srD  (l)2 _ I 

<  { 4((1  +  s2)  (<&  -  VA  -  4ft  j  +  N j)  -  4ft2]  -  {DW  +  sirDSl  +  4ft)]2 
The  above  bound  will  be  meaningless  for  those  profiles  satisfying  NA  -  0.  A  bound  for  those  cases  can 

also  be  obtained  from  Eq.  (33)  and  the  growth  rate  is  then  bounded  by 


,  _  f  IDW  +  sirDd  +  4ft)]2  _  ... 

“?<| - 4(TT7) - 


In  the  following  we  will  further  investigate  upper  bounds  for  possible  unstable  waves  and  compare 


them  with  the  stability  characteristics  previously  obtained.  Three  special  cases  will  be  discussed. 


A.  The  axisymmetric  case  (m  «■  0) 

For  mathematical  simplicity  in  this  case,  we  further  ignore  the  axial  magnetic  field,  i.e.,  WA  -  0. 
The  integral  equation  for  the  present  case,  reduced  from  Eq.  (28),  becomes 

f  p0(\D*4>\2  +  k2\+\2)rdr  +  J pok2  |$  -  VA  -  rdr  -  0  (60) 

as  the  governing  relation  for  possible  instabilities.  From  Eqs.  (20),  (25)  and  (60)  we  follow  the  pro¬ 
cedures  used  by  Kochar  8c  Jain  (1979)  in  their  derivation  of  the  semi-ellipse  theorem  in  two- 
dimensional  stratified  flows  to  obtain 


where  Jk  as  given  in  Eq.  (39)  is  restricted  to  be  less  than  1/4  if  Eq.  (61)  is  valid.  The  subscripts  max 
and  min  represent  the  maximum  and  minimum  of  the  quantities  within  the  flow  domain.  The  semi¬ 
ellipse  theorem  for  the  axisymmetric  case  thus  states  that  the  unstable  axial  phase  velocity  must  lie 


within  a  semi-ellipse  in  the  complex  phase  velocity  plane  as  described  by  Eq.  (61).  The  semi-ellipse  is 
bounded  by  the  upper  and  lower  limits  of  the  axial  flow  velocity  and  is  exactly  the  same  the  one  in 
two-dimensional  stratified  flow  except  that  the  local  Richardson  number  is  defined  as  in  Eq.  (39) 
instead. 

B.  The  azimuthal  case  ( k  «  0) 

By  using  the  Boussinesq  approximation,  Acheson  (1973)  was  able  to  demonstrate  the  westward 
drift  phenomenon  in  his  study  on  hydromagnetic  wavelike  instabilities  in  a  rapidly  rotating  stratified 
fluid.  In  addition,  a  quadrant  theorem  reminiscent  of  the  semi-circle  theorem  in  two-dimensional 
stratified  flows  was  also  derived  for  slow  amplifying  waves.  Even  though  his  criteria  were  obtained  for 
nonaxisymmetric  modes  under  certain  assumptions,  any  attempt  to  infer  that  the  azimuthal  modes  carry 
the  same  characteristics  as  the  nonaxiymmetric  ones  is,  at  best,  uncertain.  In  the  present  case,  we  will 
investigate  this  uncertainty.  The  procedure  to  be  used  to  construct  the  bound  for  instabilities  is  similar 
to  the  one  used  by  Fung  (1982)  in  his  study  of  nonaxisymmetric  instability  on  rotating  flows  except 
that  the  present  case  is  complicated  by  the  magnetic  field. 


From  the  derivation  of  the  sufficiency  condition  for  stability,  the  integral  equation  of  unstable 
waves  for  k  —  0  reduced  from  Eq.  (28)  is  written  as 


J*  po(rJ|Z)*<6l2  +  m2\<t>\2)rdr 
+  j  po«2^n  5 


Itfl2 

f  Pom2 

J  IaM2 


4-t 


D\+\- 
1 


l  mrDCKmSl  —  w.) 
7+  pTp 


rdr 


1  + 


m2  nj 
I*!2 


r2(Z>a)2|<6l2rdr  —  0. 


Substituting  Eq.  (25)  into  (62)  and  applying  the  Schwarz  inequality  lead  to  the  inequality 
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+ 1 

V  Cmi 

Next  we  separate  the  real  and  imaginary  parts  of  Eq.  (50)  for  instability  into 

f  (n2  -  2  cmtl  -  ni)  Y  mdr  +  (d  -  4)J  Xmdr  -  0 


(62) 


(63) 


(64) 
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and 


f  (lYmdr-cm,fxmdr-0.  (65) 

Let  Omx  be  the  upper  bound  and  ftmin  be  the  lower  bound  of  the  angular  velocity  within  the  flow 

held.  Because 

/  (ft  -  (W  (ft  -  ftmin)  Y mdr  <  0  (66) 

incorporating  Eqs.  (64)  and  (65)  into  (66)  yields 

cmi  ~  (ft  m*x  ^  ^min^  cmr  ftmix^min  (O^^min)  £  ^-m^r 

+  /  Ifim«flmin(^Po)  +  Z>(p0fti»'Wdr  <  0.  (67) 

Equation  (67)  implies  that  the  complex  angular  phase  velocity  will  no  longer  be  bounded  by  a  semi* 

circle  if  ft,„„ftmjnZ>po  +  Z)(p0ft \)  <  0.  This  characteristic  will  be  demonstrated  by  an  analytical  solu¬ 
tion  to  the  stability  equation  to  be  given  in  the  next  section.  To  construct  a  bound  for  possible  unstable 
waves,  we  assume  that 

^  m»x^  min^PO  ^  ®  (68) 

and  combine  Eq.  (63)  and  (67)  to  obtain 


A 

+  rZ>(p0ft  j) 

Po(rftfl)2 

min  | 

T 

F 

V  Cmt 

r 2 
cm4 


< 


(flJJmiir 


(69) 


Thus  the  complex  angular  phase  velocity  for  unstable  waves  will  be  bounded  by  a  curve  described  in 
Eq.  (69).  Another  conclusion  we  can  directly  draw  from  Eq.  (69)  is  that  stabilities  will  occur  if  the 
minimum  absolute  value  of  the  angular  Alfvin  velocity  exceeds  half  of  the  maximum  angular  velocity 
difference,  i.e., 

ln„U>  }  (ftmu-Om*).  (70) 
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One  should  keep  in  mind  that  both  criteria,  Eq.  (69)  and  Eq.  (70),  are  valid  only  with  the  restriction 
given  in  Eq.  (68).  The  bound  described  by  Eq.  (69)  depends  not  only  on  the  density  but  also  on  the 
upper  and  lower  bounds  of  the  angular  velocity  and  of  the  Alfv6n  waves.  The  bound  in  Eq.  (69) 
reduces  to  the  semi-ellipse  theorem  derived  by  Fung  (1982)  should  the  azimuthal  magnetic  field  be 
deleted  from  the  flow.  The  bound  in  Eq.  (69)  also  indicates  that  the  amplifying  waves  do  not  neces¬ 
sarily  propagate  against  the  basic  rotation.  Accordingly,  the  westward  drift  that  exists  in  the  arbitrary 
amplifying  modes  does  not  prevail  in  the  azimuthal  case. 

C.  Uniform  rotation  and  constant  angular  Alfven  wave 


A  solution  for  azimuthal  modes  described  by  a  semi-circle  can  be  derived  for  all  unstable  waves  of 
the  flow  with  uniform  rotation  and  constant  angular  Alfven  velocity  regardless  of  the  particular  form  of 
the  density  distribution.  For  ft  -  fl0  and  —  ftA0,  Eq.  (46)  under  the  present  assumption  yields  a 
simple  form 


where 


D(p^r2D*u)  —  (Ar(Z)po)  +  m2p<>iu  —  0 


(71) 


and 


A  *  Ar  +  i\t 


(n  0-O2-n20 


A'  “  77S - b ?rriT  N(n0  -  C* r)2  -  -  <£]  KOo  -  -  °io  -  4J 

l(n0-  cj2-  n^or 

+  4n0  (ftp  -  <w)  ciJ 


(72) 


A, 


_ _ 

Kfto-  “  O^ol2 


flip 

n<? 


Equation  (71)  and  the  conditions  that  u  vanishes  at  the  inner  and  outer  boundaries  form  a  Sturm- 
Louiville  system  having  the  following  characteristics:  (1)  A  is  always  real  indicating  A(  —  0,  and  (2)  A 
and  Dpo  are  of  opposite  signs.  These  characteristic  can  also  be  shown  by  applying  the  integral  method 
to  Eq.  (71)  or  directly  obtained  from  Eq.  (50)  under  the  present  assumptions.  The  first  characteristic 
states  that  the  complex  angular  phase  velocity  for  all  unstable  waves  must  lie  on  a  semi-circle  described 


by 
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and  since 


A  “  A,  “ 


(ni-  nio) 


i-tF-  >o 


(74) 


l(n0-  O2  -  n^ol2  l  ftoj 

must  hold  for  all  unstable  waves,  the  second  characterisitic  clearly  demonstrates  that  instabilities  are 
impossible  when  Dp0  >0.  As  a  matter  of  fact  this  conclusion  can  also  be  drawn  directly  from  Eq.  ($1) 
even  for  arbitrary  angular  Alfven  velocity.  Instabilities  of  this  type  are  certainly  of  centrifugal  origin 
since  no  shear  layers  exist  within  the  flow  held.  In  the  absence  of  the  azimuthal  magnetic  force,  the 
present  semi-circle  for  unstable  waves  reduces  to  the  one  obtained  by  Fung  (1982)  in  his  analysis  on 
nonaxisymmetric  instability  of  vortex  flows.  Once  again,  the  unstable  waves  do  not  propagate  against 
the  basic  rotation  unless  >  Oq- 


SOME  EXACT  SOLUTIONS 


To  demonstrate  the  validity  of  the  sufficiency  condition  for  stability  and  of  the  bounds  for 
unstable  waves,  we  proceed  to  construct  several  flow  profiles  and  obtain  exact  solutions  to  the  govern¬ 
ing  stability  equation.  Three  types  of  vortex  flows  will  be  considered. 

(I)  Uniform  Rotation 


First  consider  a  flow  profile  with  the  distribution  as  follows: 


n(r)-n2  n<4(r)-n(42 

JF(r)-»y2  fVA(r)  -  WA2.  (75) 

p0(r)  “  pi(r/R)* 

Here  R,  <r  and  all  the  quantities  with  the  subscript  2  are  arbitrary  constants.  The  perturbation  velocity 
in  the  radial  direction  obtained  by  solving  Eq.  (20)  is 


..  -t-'I.U  .  2».w!n,-tf,1n,,> 

— 

.  »[<r  .  2m(W2n2  -  NA2ClAi) 

+  *  T+  zrrw, 


kgr/',(kgr) 

l,(kgr) 


/,(**r) 


ktrK’Mar ) 
Kt(kgr) 


K,{kgr) 


(76) 


with  the  corresponding  perturbation  pressure  equal  to 


where 


Q  2 


-Hi*}  - 


9  y 

—  r  2  lAIy(kgr)  +  BK,{kgr)\ 

A 


(77) 


,  <r(n22-  n^2)  . 

q  -  1 - — - : - 4 


S\  -  N2a1 


^V2fl2  —  NA2ClA2 


Ni  -  N}2 


-Uv  + 


f+ 


2m(A(2ft2  — 


Atf  -  N}2 


1  J/2 


iV^  kW j  4*  ntd}  —  w 

NA2  "  ktVA2  +  A2 

and  Ir(kgr)  and  K,(kgr)  are  the  modified  Bessel  functions  of  the  first  and  second  kind  of  order  v.  The 
prime  denotes  the  total  derivative  with  respect  to  the  arguments  shown.  We  further  consider  the 
present  flow  to  be  confined  within  two  solid  boundaries  located  at  R\  and  R2.  The  secular  relation 
obtained  by  applying  the  boundary  conditions  u2(R\)  —  u2(R2)  -  0  is 
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(78) 


where 


*]  -  kRj  j -  1,  2. 

It  is  difficult  to  observe  the  general  behavior  of  Eq.  (78)  because  both  the  argument  and  the  order  of 
the  modified  Bessel  functions  involve  the  complex  eigenfrequency.  However,  stability  behaviors  for 
some  special  cases  can  still  be  observed  and  compared  with  the  earlier  obtained  criteria.  One  such  case 
is  the  axisymmetric  mode  with  zero  axial  magnetic  field.  The  asymptotic  expansion  of  the  BesMl  func¬ 
tions  for  large  axial  wave  numbers  allows  us  to  obtain  the  explicit  solution 


19 


-  Br  *  (<r  +  4)Oj-<rOj2  V* 

Ck  "  k  "  2  A*  +  aV/(Aj  -  *,)J 


where 


a  -  0.  ±  1.  ±2 . 

Equation  (79)  shows  the  Sow  will  be  stable  if  and  only  if 

(o-  +  4)nf-<rnJ2  >  0.  (80) 

This  is  exactly  what  the  generalized  Michael  condition  predicts  in  Eq.  (40).  The  stability  domain  is 
plotted  in  Fig.  1.  For  axisymmetric  modes,  the  presence  of  the  azimuthal  magnetic  field  stabilizes  the 
flow  for  negative  density  gradients  and  destabilizes  the  flow  for  positive  density  gradients. 

Another  special  case  is  the  one  for  azimuthal  modes.  The  solution  to  Eq.  (78)  for  k  —  0  is  found 


(<r  +  A)fl2  ±  -\/(tr  +  h)  (trflfi  +  Afi^) 


where 


^  +|in(K!/Jl2)  j 

and  a  is  an  arbitrary  integer.  Because  (<r  +  A)  is  always  positive  definite,  the  flow  will  be  stable  if 


0’°i  +  |"H  +  T-+  ln(J?,/JtI)j  }°j  »  °' 


u  p  -  4  |  me*,/*,)  |  ' 

Unlike  the  stability  behavior  in  the  axisymmetric  case  the  angular  Alfvin  velocity  always  stabilizes  the 
azimuthal  disturbances.  This  conclusion  is  consistent  with  the  stability  characteristics  predicted  by  Eq. 
(51)  for  azimuthal  modes.  The  stability  boundaries  described  in  Eq.  (82)  are  plotted  in  Fig.  2.  Due  to 
the  presence  of  the  azimuthal  magnetic  field,  the  flow  can  still  be  stable  even  when  the  density  is 
decreasing  radially  outwards.  Another  observable  stability  characteristic  in  the  case  of  azimuthal  modes 
is  that  the  solution  for  unstable  waves  in  Eq.  (81)  ties  on  the  semi-circle  as  prescribed  by  Eq.  (73). 


(2)  Tw*-Rcgi«i  Flew 

The  second  flow  profile  we  would  like  to  investigate  is  m  two-region  flow  with  flu  iatsrfftps 
between  the  inner  and  miter  region  located  at  r  •  it.  The  flow  distributions  are  as  follows: 


LVJttS&v  ft*  Tv? 


m 


Fig.  I  —  The  stability  domain  for  axisymmetric  modes 


Fig.  2  —  The  stability  for  azimuthal  modes 
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ft(r)-  ft, 

W(r)  -  Wj 
fl^(r)  “ 
f^(r)-  WA 
p0(r)  -  Pj 

where  all  the  quantities  with  index  j  are  constants.  After  matching  the  interfacial  conditions  in  Eqs. 
(23)  and  (24),  the  secular  relation  for  stability  is  found  to  be: 
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If  the  axial  Alfven  velocities  are  neglected  and  the  perturbations  are  constraavwa  to  be  axisyi&metric, 
the  asymptotic  expansion  of  the  modified  Bessel  functions  for  large  axial  wave  numbers  allows  us  to 
solve  Eq.  (83)  and  the  resultant  condition  for  stability  is 


^-  +  ^■1  !p2(ni-n2I)-p,(fl?-n3i)]-*2(»'l-  w^o. 


(84) 


This  condition,  composed  of  the  centrifugal  force  jump  and  the  axial  velocity  difference  at  the  interface, 
is  consistent  with  Eq.  (43)  under  the  same  assumption.  The  flow  will  be  unstable  if  the  centrifugal 
force  at  the  interface  is  decreasing  outwards.  The  axial  velocity  difference  always  destabilizes  the  flow 
except  for  Wx  -  W2.  In  that  case,  the  generalized  Michael  condition  in  (40)  is  recovered  as 


p2(ftf-  ftj2)  -Pi(ftj*-  ftj,)  >  0 


(85) 
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for  stability.  As  a  matter  of  fact,  condition  (85)  can  easily  be  obtained  by  integrating  Eq.  (40)  across 
the  interface,  since  the  flow  is  stable  within  the  inner  and  the  outer  regions. 


Another  solution  we  would  like  to  obtain  from  Eq.  (83)  is  one  for  the  azimuthal  modes.  For 
k  -  0,  Eq.  (83)  yields  a  simple  solution 


where 


(m  —  Dpifli  4-  (m  4-  l)p2fl2  ±  V0 

m(pi  +  p2) 


(86) 


0  —  — (m2  —  l)p,p2(n,  -  fl2)J  +  (p2  -  pi)  Km  -  Dpifl2  +  (m  +  l)p2flf] 

+  m(pi  +  p2)  Km  —  Dpifl^i  +  (m  +  l)p2flj2I.  (87) 

Stability  of  the  flow  will  be  guaranteed  if 

0  >  0.  (88) 
As  in  our  previous  discussion  of  Eq.  (51)  for  azimuthal  modes,  the  first  term  in  Eq.  (87)  is  the  shear 

influence  generated  by  the  difference  in  the  angular  velocities  at  the  interface  which  always  destabilizes 

the  flow.  The  second  term  represents  the  effect  of  the  density  variation  in  the  centrifugal  force  field 

created  by  the  rotation  of  the  fluids  in  the  inner  and  outer  regions.  Stabilization  effect  requires  the 

density  gradient  at  the  interface  to  be  positive.  The  last  term  in  (87)  is  the  effect  of  the  azimuthal 

magnetic  field  which  always  stabilize"  the  flow. 


For  instability,  a  semi-circle  bound,  if  valid  for  the  present  profile,  would  read 


Cm.  - 


n,  +  n2 


+  'ml  < 


n2  —  n  i 


(89) 


Substituting  Eq.  (86)  into  Eq.  (89),  we  find  that  the  above  inequality  will  be  satisfied  if 


(p2  -  p,)n,n2  +  Km  +  l)p2Cl  22  +  (m  -  Dpifl  j(]  >  0.  (90) 

This  satisfies  the  testriction  in  Eq.  (68)  for  the  validity  of  a  semi-circle  bound.  Equation  (90)  also  indi¬ 
cates  that  unstable  waves  will  no  longer  be  bounded  by  a  semi-circle  if  the  inequality  is  violated. 
Unlike  the  semi-circle  bound  on  all  unstable  waves  in  two  dimensional  stratified  flows,  the  semi-circle 
theorem  in  vortex  flows  does  not  in  general  provide  a  bound  on  unstable  waves  if  the  restriction  in  Eq. 
(68)  is  violated. 
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(3)  Three-Region  Flow 


As  the  last  example  to  demonstrate  the  criteria  previously  derived  for  uniformly  rotating  flows 
subject  to  azimuthal  perturbations,  we  consider  a  three-region  flow  with  constant  angular  velocity  fl0, 
constant  angular  Alfvln  velocity  and  a  density  distribution  specified  by 

,  0  <  r  <  /?, 

/?,  s?  r  <  R2  . 

/?2  <  r  <  <*> 

Here  p,,  p2,  p3  and  <r  are  constants  that  characterize  the  density  profile  in  the  three  flow  regions  with 
their  common  interfaces  at  the  radial  positions  R\  and  R2-  The  axial  components  of  the  velocity  and 
magnetic  field  in  this  case  can  be  arbitrary.  With  the  solutions  for  the  velocity  and  pressure  perturba¬ 
tions  given  in  Eqs.  (76)  and  (77),  we  obtain,  after  using  the  matching  conditions  in  (23)  and  (24),  the 
secular  relation  for  stability  as 

mp\  +  ^TPi  +  (P2~  Pi)A  -  p2v  mp3-  -f-p*2  +  (p3-  P*)A  +  p*2v 

- £ - R2v  - £ - R2 *  (91) 

mp\  +  yP2  +  <P2  -Pi)A  +p2v  wp3  -  y p\  +  (p3  -P*)A  -p\v 

where  A  is  given  in  Eq.  (72),  v  is  given  in  Eq.  (77),  and  p  *  *  p2(/?2/R  i)'  is  the  density  of  the  middle 
region  evaluated  at  r  -  R2.  Two  special  cases  for  the  present  flow  profile  will  be  investigated. 


Pi 

Po(r)  -  p2 

*i 

Pi 


The  first  special  cases  is  for  <r  ~  0,  i.e.,  the  fluid  is  homogeneous  in  the  middle  region.  Equation 


(91)  yields 


A  1  -  —  1  -  —  A2- m(l  +  A)  1  -  —  A  +  m2  l  +  —  +S—  +  —  -0  (92) 

Pi  Pi  Pi  |  Pi  Pi  Pi 


where 


Rjm  -  R}m 


R}m  +  R}m' 

There  are  four  general  types  of  density  profiles  for  the  flow  under  consideration.  They  are:  (a) 
Pj>  Pi  >  Pi,  (b)  pi<  pi<  pi,  (c)  pj  <  Pi  >  pi,  and  (d)  pj  >  p2  <  pi-  Combining  Eqs.  (72)  and 
(92),  we  find  that  instabilities  will  occur  when 
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and 


O^lyl  +  ft^oOxI  ±  Vx*  -  If  |)  <  O  for  type  (a) 
OjM  -  ft^oOxI  ±  Vx5  -  If |)  >  0  for  type  (b) 


(93) 


~  V7  +  W  >  0  for  types  (c)  and  (d) 
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Since  the  sums  of  the  terms  inside  the  radicals  in  Eq.  (93)  are  always  positive;  we  can  immediately  con¬ 
clude  from  the  conditions  in  Eq.  (93)  that  the  presence  of  the  azimuthal  magnetic  held,  as  in  our  previ¬ 
ous  discussion  of  Eq.  (SI),  always  stabilizes  the  flow.  It  is  obvious  that  the  instability  condition  for 
type  (a)  can  never  be  satisfied,  as  also  predicted  by  Eq.  (51).  since  the  density  is  increasing  radially  out¬ 
wards.  In  the  absence  of  the  azmuthal  magnetic  field,  the  conditions  in  Eq.  (93)  reduce  to  those  in 
Fung  (1982)  and  the  flow  is  always  unstable  except  for  the  density  profile  in  type  (a).  It  can  also  be 
shown  that  for  instability,  all  the  unstable  waves  for  the  density  profile  in  types  (b),  (c)  and  (d)  lie  on 
the  semi-circle  described  by  Eq.  (73). 


The  second  special  case  we  would  like  to  examine  is  a  continuously  varying  density  distribution. 
Equation  (91)  forp)  -  p2  and  pj  —  pi(RjRi)9  yields  the  simple  form 

(2m2  +  o-A  —  2mv)Rfy  —  (2m2  +  <rA  +  2  mv)R}r.  (94) 

Equation  (94)  has  the  same  form  as  the  one  (Eq.  22)  in  Fung  &  Kurzweg  (197S)  except  for  the  terms 

involving  A  and  v  due  to  the  presence  of  the  magnetic  field.  A  simple  solution  to  Eq.  (94)  for  v  —  0 

allows  us  to  solve  for  the  complex  phase  velocity  as 


Since  the  sum  of  the  terms  in  the  first  bracket  inside  the  radical  is  positive,  we  can  again  conclude  the 
azimuthal  Aflven  velocity  also  stabilizes  the  flow.  Furthermore,  the  complex  phase  velocity  for 
unstable  waves  in  this  case  also  lies  on  a  semi-circle  as  predicted  by  Eq.  (73). 

CONCLUSIONS 

Some  general  stability  criteria  for  a  general  type  of  vortex  flows  of  conducting  fluids  with  axial 
velocity  components  under  the  influence  magnetic  field  are  derived.  Exact  solutions  to  the  governing 
stability  equation  for  some  special  flow  profiles  are  obtained  and  compared  with  the  earlier  derived  cri¬ 
teria  of  a  less  general  nature. 

The  derived  sufficiency  condition  for  stability  is  generally  unseparable  from  the  complex  eigenfre- 
quency  because  of  the  presence  of  the  axial  and  the  azimuthal  magnetic  field.  To  further  investigate 
the  roles  played  by  the  density,  velocity  and  magnetic  field  in  flow  stability,  three  types  of  perturbation 
conditions  are  investigated  and  compared  with  the  sufficiency  condition.  It  was  shown  that  densities 
that  increase  radially  outwards  always  have  a  stabilizing  effect.  Unlike  the  axial  velocity  which  always 
destabilizes  the  flow,  the  angular  velocity  plays  a  dual  role  in  flow  stability.  While  the  gradient  of  the 
velocity  generates  shear  effects  which  destabilize  the  flow,  the  rotation  of  the  velocity  creates  a  centrifu¬ 
gal  force  field  which  stabilizes  or  destabilizes  the  flow  depending  on  the  sign  of  the  density  gradient.  If 
perturbations  to  the  flow  are  restricted  to  be  azimuthal,  the  magnetic  field,  regardless  of  its  detailed  dis¬ 
tribution,  always  stabilizes  the  flow. 

Several  bounds  on  unstable  waves  are  also  obtained  and  compared  with  some  exact  solutions  to 
the  stability  equation.  For  axisymmetric  instabilities,  a  semi-ellipse  theorem  is  proved  to  be  valid  in  the 
absence  of  the  axial  magnetic  field.  For  azimuthal  instabilities,  an  upper  bound  on  the  complex  phase 
velocity,  reminiscent  of  the  semi-ellipse  theorem  in  the  non-magnetic  case,  is  derived  with  a  restriction. 
Such  a  restriction  indicates,  as  supported  by  an  exact  solution  to  the  stability  equation,  that  not  all 
unstable  waves  are  bounded  by  a  semi-circle  in  the  complex  phase  velocity  plane.  For  flows  with  uni¬ 
form  rotations  and  constant  Alfvin  velocities  subject  to  azimuthal  perturbations,  instabilities  can  only 


occur  when  the  density  decreases  radially  outwards.  Furthermore,  the  phase  velocity  for  such  instabili- 
ties,  regardless  of  the  detailed  density  distribution,  must  lie  on  a  semi-circle  in  the  complex  velocity 
plane.  The  semi-ellipse  and  semi-circle  bounds  derived  for  the  azimuthal  modes  clearly  show  that  the 
amplifying  waves  do  not  necessarily  propagate  against  the  basic  rotation.  As  a  conclusion,  therefore, 
the  westward  drift  derived  for  arbitrary  nonaxisymmetric  modes  does  not  generally  prevail  in  the 
azimuthal  case. 
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